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The  Department  of  Engineering  Science  at  Oxford  University,  Oxford, 
U.K.,  has  had  an  active  turbomachinery  research  group  for  over  15 
years.  Major  contributions  to  the  literature  of  fundamental 
measurements  in  turbine  blade  heat  transfer  are  well  established 
worldwide.  Professor  UaGraff  has  closely  followed  the  contributions  of 
the  Oxford  group  for  many  years,  talcing  a  sabbatical  leave  at  Oxford 
from  Syracuse  University  for  seven  months  in  1983.  The  research 
conducted  during  this  visit  lead  directly  to  a  proposal  for  a  follow-on 
experimental  program  ,  viz.,  examining  in  some  detail  the  transitional 
behavior  of  the  boundary  layer  under  the  influence  of  wake/shock 
induced  unsteadiness.  Preparation  for  this  work  began  in  the  summer  of 
1985,  with  travel  support  from  the  National  Science  Foundation  (Grant 
No.  INT-8509407)  .  The  first  stage  of  the  work  reported  herein  was 
completed  in  the  summer  of  1986.  Prof.  LaGraff  was  on-site  from  May 
to  August  1986  when  access  to  the  primary  test  facility  was  scheduled. 
Brief  visits  before  and  after  the  primary  test  times  were  also  made  in 
October  1985,  March  1986,  September  1986,  and  October  1986. 
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X .  Introduction 


A  new  type  of  transient  cascade  has  been  developed  at  Oxford 
University  which  has  substantially  improved  knowledge  of  the  distribution 
of  heat  transfer  rates  over  turbine  bladings*  This  facility,  described 
briefly  below,  allows  the  testing  of  actual  blade  profiles  at  full  scale 
engine  Reynolds  numbers,  Mach  numbers  and  gas-to-wall  temperature  ratios* 
The  theory  of  operation  and  typical  results  are  discussed  in  detail  in 
References  1  and  2.  A  summary  of  the  facility  and  instrumentation  are 
included  in  the  next  two  parts  of  this  report  for  completeness. 

(See  e.g.  Fig.  1) . 

Of  special  interest  in  comparing  cascade  tests  with  real  engine  rigs 
is  to  account  for  the  flow  unsteadiness  associated  with  wake-passing 
phenomena  between  adjacent  stages.  Some  recent  work  has  been  reported 
which  has  begun  to  unravel  this  complex  flow  field  in  a  systematic  manner 
(see  e.g.  Refs.  3  and  4).  Very  recently  Doorly  (Ref.  5)  has  reported 
detailed  time-resolved  heat  transfer  and  aerodynamic  measurements  of  a 
wake  passing  experiment  in  the  Oxford  cascade.  In  this  experiment  the 
wakes  were  generated  by  a  rotating  set  of  cylindrical  wires  upstream  of 
the  stationary  rotor  cascade  (Fig. 2) .  Significant  changes  in  the  cascade 
flow  and  heat  transfer  rates  were  observed  during  each  wake  passing  event . 
Subsequent  work  in  the  Oxford  cascade  on  a  new  blade  profile  was 
systematically  undertaken  at  various  Mach  numbers,  Reynolds  numbers  and 
blade  incident  angles  to  determine  the  wake  passing  vs.  no  wake  results 
(Ref.  6)  with  the  intention  of  ultimately  comparing  this  with  tests  on  the 
same  blade  on  the  new  MIT  Gas  Turbine  Laboratory  transient,  rotating, 
blow-  down  rig.  These  comparisons  are  not  complete  at  this  date. 
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II.  Thm  Experimental  Technique 

«) Ifltntropic  Light  Platon  Tunnel  Cascade 

Transient  techniques  have  long  been  recognized  as  a  powerful  technique 
for  heat  transfer  measurements .  (Ref.  7.).  Short  duration  shock  heated 
facilities  had  particular  success  simulating  the  high  enthalpy 
environment  characteristic  of  hypersonic  reentry  conditions  from  orbital 
and  suborbital  flight.  Inferring  instantaneous  heat  transfer  rates  from 
the  transient  response  of  the  surface  temperature  of  models  (see  next 
section)  greatly  eased  the  instrumentation  problems  in  the  high 
temperature  environment  of  steady  state  facilities,  and  substantially 
reduced  the  operating  costs  compared  to  long  duration  testing. 

In  recent  years  (Refs.  8,9,10)  transient  facilities  have  found 
application  at  the  more  modest  enthalpy  levels,  but  experimentally 
awkward  environment,  of  gas  turbine  engines  (transonic,  rotating, 

Rt  -  10*  +,  Tq  /T^  -  1.5  -4  2.0),  A  particularly  productive  aporoach  has 
been  the  one  developed  at  Oxford  University  (Ref.  1)  where  shock  heating 
has  been  replaced  by  a  piston  driven  isentropic  compression  process.  A 
fundamental  advantage  of  this  approach  is  that  the  running  time  is 
increased  by  one  or  two  orders  of  magnitude  over  the  shock  heated  mode 
while  still  offering  the  advantages  of  transient  techniques.  The 
increased  running  time  offers  significant  instrumentation  advantages  not 
possible  in  shock  heated  facilities,  including  It  greater  ease  of 
simultaneous  direct  digital  sampling  of  multichannel  instrumentation, 

2) the  possibility  of  multichannel  surface  pressure  measurements  through 
relatively  long  lengths  of  tubing  to  transducers  mounted  external  to  the 
small  models,  3)  surface  oil  drop  flow  visualization  and,  4)  airfoil  wake 
pressure  traverses  yielding  the  loss  measurements  which  are  of 


significant  interest  to  turbine  aerodynamic ists .  Light  piston  isentropic 
compression  tunnels,  while  now  well  established  at  two  centers  in  Europe 
(Oxford  and  Von  Karman  Institute),  have  not  received  much  attention  in 
the  USA.  A  comprehensive  discussion  of  the  theory  of  the  light  piston 
facility  can  be  found  in  Reference  2.  Basically  the  running  time  of  the 
piston  facility  is  limited  by  the  driven  tube  volume  and  not  its  length 
as  in  a  shock  tube.  Thus,  short  length  but  large  diameter  tubes  can 
sustain  a  large  mass  flow  even  through  the  relatively  large  throat  areas 
of  five-blade  cascade  rows  and  not  be  limited  by  the  wave  reflection 
process  that  dominates  shock  tunnels.  As  a  result,  exceptionally  steady 
flow  rates  for  up  to  0.5  seconds  are  easily  achieved  in  the  light  piston 
facility.  The  Oxford  facility  is  routinely  operated  at  full  scale  engine 
Reynolds  number  and  Mach  number  and  at  gas  to  metal  (wall)  temperature 
ratios  of  1.5  or  higher,  simulating  the  important  parameters  for  turbine 
blade  aerodynamics  and  heat  transfer.  Equally  important  is  the  fact  that 
these  parameters  can  be  varied  independently.  A  brief  discussion  of  the 
heat  transfer  instrumentation  for  the  reported  work  follows  in  the  next 
section. 

b)  Instrumentation 

The  instrumentation  central  to  the  reported  research  involves  the 
high  bandwidth  heat  transfer  measurements.  The  sensors  are  thin  film 
resistance  thermometers  fired  onto  a  blade  profile  which  has  been 
machined  out  of  a  glass  ceramic.  The  surface  temperature  is  monitored 
and  easily  converted  to  an  instantaneous  heat  transfer  rate  using  as  a 
signal  processor  an  electrical  analogue  of  the  1-D  heat  conduction 
equation.  Figure  3  shows  an  example  of  an  inst rumentated  blade,  the 
signal  processing,  and  typical  results  fc:  heat  transfer  and  surface 
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temperature  response.  For  convenient  comparisons,  it  is  usual  to 
extrapolate  instantaneous  heat  transfer  rates  for  various  points  on  the 
blade  to  a  uniform  {isothermal)  initial  surface  temperature.  Such  an 
extrapolation  procedure  is  also  shown  in  Figure  3. 

Special  electronics  have  been  reported  by  Oldfield  et.al.  (Reference 
11)  which  raise  the  effective  bandwidth  of  the  instrumentation  package  to 
100  KHz,  sufficient  for  observing  individual  wake-passing  events.  Figure 
4  shows  typical  results  for  surface  heat  transfer  from  Reference  6.  The 
high  time  resolution  data  were  recorded  on  two  8  channel  Datalab  Model 
DL2800  Transient  Recorders  where  4096  point  were  recorded  on  each  channel 
at  a  rate  of  up  to  2  MHz/channel. 

III.  Research  Proposed  on  Boundary  Layer  Transition  with  Wake  Interactions 

Recently  studies  have  been  made  of  the  effect  of  the  unsteady  flows 
caused  by  nozzle-wake/rotor  blade  interactions.  (Refs.  5  and  6).  These 
studies  have  shown  that  several  effects  occur  in  an  interactive  manner  and 
include  boundary  layer  disturbances  due  to  incidence  changes,  shock 
impingement  (from  the  NGV  wake)  and  early  transition  caused  by  the 
increase  in  free  stream  turbulence  as  the  wake  from  the  NGV’s  passes 
through  the  turbine  passage.  In  previous  work,  tracking  of  the  boundary 
layer  transition  has  been  done  with  the  span-wise  thin  film  transducers 
approximately  10mm  wide  which  ’average*  the  heat  transfer  rate  spanwise. 

The  transducers  themselves  have  a  very  wide  bandwidth  0.1  Hz  100  kHz  and  so 
can  resolve  the  rapid  phenomena  associated  with  conventional  "turbulent 
spot"  growth,  (Refs.  12,13,14  ).  What  was  proposed  was  to  determine 
whether  the  process  of  transition  on  a  turbine  blade  resembles  that  on  a 
stationary  aerofoils.  The  technique  developed  at  Oxford  produces  a 
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relative  motion  between  a  rotor  and  nozzle  guide  vane  by  means  of  a 
rotating  set  of  bars  (20-25000  RPH)  upstream  of  the  rotor.  (Fig.  2) .  The 
wake  is  designed  to  be  correct  in  magnitude  of  momentum  deficit  and 
dimensions  to  replicate  that  of  the  NGV.  The  rotor  is  stationary.  High 
speed  heat  transfer  gauges  on  the  rotor  are  used  to  determine  the  blade- 
wake  interaction  and  are  capable  of  tracking  the  conjectured  transition 
1  spots  *  . 

The  proposed  first  year  study  was  a  basic  analysis  of  the  transition 
process  on  2-D  cascade  airfoil  under  the  influence  of  simulated  wakes 
passing  over  the  surface.  The  methodology  involved  wide  bandwidth  heat 
transfer  instrumentation  (closely  spaced  narrow  thin  film  thermometers) 
and  high  speed  digital  data  recorders  allowing  up  to  16  channels 
of  simultaneous  time-resolved  data.  With  this  data,  correlation  spectra 
could  be  extracted  while  tracking  transitional  events  in  the  blade 
boundary  layer. 

IV.  Summary  of  Results 

A  fully  documented  scientific  summary  of  the  results  of  the 
first  year  of  the  program  are  contained  in  Ref.  16  (a  copy  is  attached 
in  Appendix  A) .  This  paper  has  been  accepted  at  the  Japanese  Society 
of  Mechanical  Engineers/American  Society  of  Mechanical  Engineers  2nd 
Joint  Thermal  Engineering  Conference  scheduled  for  March  1987.  A  brief 
summary  (Ref.  17)  of  the  work  was  submitted  as  part  of  a  more 
comprehensive  paper  submitted  to  the  AGARD  Propulsion  and  Energetics 
Panel  Meeting  in  Neubiberg,  W.  Germany  in  September  1986  (A  copy  of 
this  paper  is  attached  in  Appendix  B) . 

The  results  verified  first  of  all  the  experimental  approach.  Narrow, 
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closely  spaced,  thin  film  constant  current  surface  thermometers  coupled 
with  wide  bandwidth  electrical  analogue  circuits  and  high-speed-multi* 
channel  digital  data  acquisition  could  detect  and  precisely  track 
transient  aerodynamic  events  on  the  2-D  turbine  blade  surface. 
Specifically,  discrete  turbulent  spot  events  were  seen  to  appear  in 
the  laminar  boundary  layer  heat  transfer  traces,  to  grown  and  finally 
merge  (in  the  high  Reynolds  number  case)  giving  a  "steady”  turbulent 
heat  transfer  level*  The  heat  transfer  traces  between  the  turbulent 
spots  returned  precisely  to  the  levels  of  the  laminar  values  (obtained 
over  the  entire  blade  with  low  free  stream  turbulence)  .  The  detail  of 
the  transitional  signals  allowed  individual  spots  to  be  identified  and 
followed  from  film  to  film.  Correlation  spectra  showed  very  high 
correlation  coefficients  between  channels  over  a  time  span  which 
included,  multiple  turbulent  spot  events.  Quantification  of  the  state 
of  the  boundary  layer  was  easily  achieved  using  a  standard 
intennmittency  definition. 

After  quantifying  the  state  of  the  individual  boundary  layer  {for 
both  the  laminar  case  and  natural  transition  case),  the  tests 
proceeded  to  investigate  the  influence  of  the  wake-passing  events  with 
and  without  the  shock  waves  which  would  occur  on  a  transonic  stage. 
Although  the  behavior  was  quite  complex  the  basic  results  indicated 
that  the  region  between  the  wake-passing  events  behaved  in  much  the 
same  way  as  the  natural  (undisturbed)  boundary  layer  over  the  entire 
blade.  Coincident  with  the  wake  and  shock  disturbances,  however,  were 
unsteady  excursio  s  in  the  heat  transfer  which  were  very  marked  near 
the  front  of  the  blade  (including  apparent  separation  of  the  boundary 
layer  as  the  shock  swept  by)  and  continuing  as  an  enharced  heat 
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transfer  along  the  blade  surface  as  far  as  the  instrumentation 
covered. 

Preparations  are  underway  for  the  extension  of  the  present  studies 
to  a  fully  annular  rotating  transonic  stage.  This  worlc  will  enter  the 
test  phase  in  the  summer  of  1987. 
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ABSTRACT 

Results  are  presented  illustrating  the  detailed 
behaviour  of  the  suction  surface  boundary  layer  of  a 
transonic  gas  turbine  rotor  in  a  2-D  cascade  under  the 
influence  of  both  freestream  turbulence  and  simulated 
nozzle  guide  vane  wakes  and  shocks.  The  instrumentation 
included  thin  film  resistance  thermometers  along  with 
electrical  analogues  of  the  1-D  heat  conduction 
equations  to  obtain  vide  bandwidth  heat  transfer  rate 
seasureaents  in  a  short  duration  wind  tunnel.  This 
instrumentation  provides  sufficient  time  resolution  to 
track  individual  wake  and  shock  related  events  and  also 
the  turbulent  bursts  of  a  transitional  boundary  layer. 
Wide  bandwidth  surface  pressure  transducers  and  spark 
Schlleren  photography  were  used  in  support  of  these  heat 
transfer  measurements.  The  results  showed  there  to  be  a 
direct  relationship  between  the  passage  of  wake 
disturbances  and  transient  surface  heat  transfer 
enhancements.  It  was  possible  to  track  both  wake  and 
transitional  events  along  the  surface  and  to  compare 
these  with  the  expected  convection  rates.  Analysis  of 
the  signals  allowed  direct  calculations  of  intermittency 
factors  which  compared  well  with  predictions.  Additional 
effects  due  to  a  moving  shock /boundary  layer  interaction 
were  investigated.  These  resulted  in  marked  variations 
in  beat  transfer  rate  both  above  and  below  the  laminar 
values.  Those  excursions  were  associated  with  separation 
and  re-attachment  phenomena. 

NOMENCLATURE 

a  -  Sound  speed 

*  Tangential  (or  true)  chord 

k  ■  Thermal  conductivity  (W/mK) 

4  -  Heat  transfer  rate  (W/ra*) 

«  dU 

%  -  Acceleration  parameter, 

M  -  Isentroplc  Mach  number 

Nu  *  Nusselt  number 

R  -  Radius  of  curvature  (am)  or 

the  gas  constant  <J/kgK) 

Re  -  Reynolds  number  (based  on  Inlet  total 

conditions,  isentroplc  exit  Mach  number  and 
tangential  chord) 


3  *  Blade  surface  perimeter 

t  ■  Time 

T  *  Temperature  (K) 

Tu  »  Turbulence  level  *  u*/U 

u*  «  Fluctuating  velocity  (m/s) 

U  =  Velocity  (0/3) 

V  *  Voltage 

P  =»  Pressure 

x  -  Surface  distance  from  the  leading  edge 
p  -  Gas  angle  (measured  from  the  axial  direction) 
y  =  Ratio  of  specific  heats 

y  «  Intermittency  factor 

t i  *  Viscosity  (kg/ms) 

u  »  Kinematic  viscosity  *  p/p 

p  Density  (kg/oM 

Subscripts 

-  =  Freestream 

o  *•  Total 

1  *  Inlet 

2  -  Outlet 

m  Parameter  at  the  selected  measuring  instant 

rel  *  Condition  relative  to  the  rotating  bar 

INTRODUCTION 

The  need  to  understand  the  processes  by  whief 
transition  takes  place  in  a  boundary  layer  and  it; 
susceptibility  t  *  disturbance  by  events  occurring  in  th» 
freestream  have  led  to  a  growing  amount  of  research 
much  of  it  aimed  towards  the  gas  turbine  industry.  Ar 
increase  -in  lbs  level  of  freestream  turbulence  has  beet 
shown  to  encourage  the  transition  process  (Blair,  1982 
Rued  and  vmtig,  1084)  following  the  turbulent  apo* 
models  proposed  by  Dchubauer  and  Klebanoff  (1956).  Thi 
effects  of  intermittent  changes  in  flew  conditions  sue) 
a:,  those  caused  by  the  Nozzle  Guide  Vane  (NGV )  and  roto: 
L.ade  lr.ttM  actions  have  also  been  shown  to  hav< 
significant  influences  on  heat  transfer  rate  to  th< 
blades,  as  detailed  by  Dunn  (1985)  and  Doorly  (1983) 
The  tests,  reported  In  this  paper  Investigate  both  th» 
transitional  phenomena  and  the  influent  of  a  simulate, 
NGV  on  a  rotci  Made  suction  surface  boundary  layer.  Th* 
use  of  wide  tandwidth  heat  transfer  lnstruoentatior 
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Figure  1  -  Instrumented  Model  Details  ana  Test 
Condi tions. 


enables  individual  events  such  as  turbulent  spots 
developing  in  the  transitional  boundary  layer  to  be 
easily  identified,  at  freestream  conditions 
corresponding  to  the  gas  turbine  environment  In  terms  of 
Mach  number,  Reynolds  number  and  gas-to-wall  temperature 
ratio. 

Experimental  Approach 

The  present  study  was  carried  out  on  a  transonic 
gas  turbine  rotor  blade  in  the  Isentropic  Light  Piston 
Tunnel  (XLPT).  as  described  by  Schultz  et  *1.(1977),  The 
ILPT  is  a  short  duration  facility  which  for  a  typical 
run  time  of  about  0.5  seconds  gives  a  gas-to-wall 
temperature  ratio  of  1*5  .  For  these  tests  a  5-bladed 
cascade  of  50  mm  inlet  span  was  used  with  a  Mach  number 
of  0.38  at  inlet  and  1.18  at  the  exit  plane  with  a 
corresponding  Reynolds  number  of  0.52  x  10*.  The 
freestream  turbulence  level  could  be  varied  from  a 
background  level  of  leas  than  0.8%  to  a  "high1*  level  of 
about  3%  by  use  of  a  grid  of  2  mm  bars  placed  208  am 
upstream  of  the  cascade  inlet  plane.  The  operating 
conditions  are  detailed  In  Figure  1,  using  the 
definitions  ; 

po  c 

Re  -  Reynolds  number  *  _ LI 

9 

with  viscosity  derived  from  Sutherland’s  law. 

M  •  Isentropic  Mach  number  *  U/at 


where  sound  speed  at 


-Air: 


size  smaller  than  the  film  coverage.  Heat  transfer  rate 
was  obtained  using  a  well-established  technique  (Schultz 
and  Jones,  1973)  by  usin*  electrical  analogues  of  the 
aomi -inf inite  conduction  equation  to  give  a  direct 
ttcasur^ncnt  of  heat  transfer  rate  with  a  bandwidth  of 
100  kHz,  as  described  by  Oldfield  et  al,<198<).  Hide 
bandwidth  surface-mounted  pressure  transducers  (Kulite 
type  XC0-062-50D)  were  used,  together  with  spark 
Schlieren  photographs  in  the  analysis  of  the  data. 

In  order  to  further  characterise  the  freestream 
conditions,  constant  temperature  hot  wire  anemometer 
measurements  were  made  at  the  channel  centra  line  at  the 
entrance  plane  of  the  cascade  leading  edge.  The  probe 
used  was  a  DISA  Pll  5  micron  platinum  coated  tungsten 
wire  of  2  mm  length  and  connected  to  a  DISA  55001 
constant  temperature  anemometer  bridge  tuned  to  a 
bandwidth  of  greater  than  50  kHz. 
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The  position  of  the  heat  transfer  instrumentation 
is  also  shown  in  Figure  1,  with  the  3  7  thin  film  gauges 
spaced  around  the  suction  surface.  For  these  tests  the 
gauges  were  spaced  2.5  mm  apart  freo  9%  to  50^  of  the 
surface  perimeter  of  64.2  mm,  with  a  5  mm  spacing 
thereafter.  The  gauges  were  0.5  am  wide  and  covered  4 
mm  of  the  span,  compared  with  the  10  mm  gauges  used 
previously  in  the  cascade  to  enable  smaller  three* 
dimensional  events  to  be  observed.  The  gauge 
effectively  averages  the  heat  transfer  rate  along  its 
length,  thus  reducing  the  magnitude  of  any  events  of 


Figure  2  -  The  Cascade  with  Rotating  Bar  Configuration. 
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The  effects  of  the  NGV-rotor  interaction  were 
simulated  by  2  mm  diameter  bars  attached  to  a  disk  spun 
upstream  of  the  cascade,  as  described  by  Doorly  and 
Oldfield  (1985b).  and  Illustrated  in  Figure  2.  By  this 
means  It  is  possible  to  simulate  the  inlet  perturbation 
In  the  rotor  relative  frame  for  the  cascade,  to  give  the 
correct  NGV  wake  orientation  and  scale.  This  method  was 
used  by  Ashworth  et  al.U9$5)  for  the  same  profile  at  an 
Inlet  angle  or  58*  with  a  bar-to-cascade  mid-span  pitch 
ratio  of  1.7  representative  of  a  realistic  turbine 
design.  The  present  study  reduced  the  number  of  bars  to 
2  giving  a  pitch  ratio  of  13.4  so  that  individual 
periodic  disturbances  could  be  identified,  giving  a  bar¬ 
passing  frequency  of  740  Hz.  In  the  fast  sampling 
period  of  6  milli-seconds  it  was  possible  to  capture  4.4 
bar-passing  cycles,  so  that  the  repeatability  of  events 
could  be  determined. 

The  data  was  acquired  In  several  ways: 

(i)  Slow  samples  were  taken  at  500  Hz  using  a  64- 
channel  DataTranslatlon  DT1711  A-D  board,  with  a 
PDP11/34  controlling  this  acquisition  and  the 
tunnel  housekeeping  channels  establishing  the 
steady  flow  conditions.  From  this  data  mean 

values  of  heat  transfer  rate  are  deduced  by 
choosing  a  time  interval  around  a  measure  point  in 
the  run  for  which  conditions  are  steady  (Oldfield 
et  al..  1978). 

(11)  A  Datalab  DL2800  transient  recorder  was  used  to 
acquire  16  channels  of  4096  points  at  a  sampling 
rate  of  500  kHz  at  a  time  in  the  run  around  the 
measure  point.  These  results  form  the  basis  of 
the  transient  analysis  in  this  paper. 

(ill)  Schlleren  photographs  stored  a  "snapshot"  of  the 
flowfield  as  an  aid  in  interpretation  of  the 
surface  measurements. 


STEADY  BASELINE  FLOW  CONDITIONS 

Measurements  of  the  mean  aerodynamic  and  heat 
transfer  properties  of  the  rotor  cascade  are  given  in 
Figure  3  for  the  suction  surface  at  the  test  condition 
of  the  present  experiment.  The  he^t  transfer 
coefficient  given  as  Nusselt  number  is  seen  to  increase 
at  and  n«.ur  the  leading  edge  as  expected  with  higher 
freestream  turbulence  and  blade  Reynolds  numbers  (Lowery 
and  Vachon.  1975).  As  the  boundary  layer  develops  along 
the  surface  the  heat  transfer  rate  in  the  high  and  low 
freestream  turbulence  cases  merge  together  at 
essentially  the  laminar  level.  The  low  turbulence  case 
heat  transfer  rate  (Tu  <  0.8%)  remains  laminar  to  the 
x/s  position  of  0.8  whereas  in  the  high  turbulence  case 
(Tu  *  3%)  begins  to  Inorease  in  proportion  (as  will  be 
shown  in  the  following  section)  to  Increasing  levels  of 
turbulent  spot  intermi ttency .  In  the  high  Reynolds 
number  (and  high  freestream  turbulence)  case  the  heat 
transfer  rate  never  comes  down  to  the  laminar  level  and 
begins  climbing  due  to  transitional  intermittency  soon 
after  the  x/s  *  0.15  location.  The  freestream  velocity 
as  determined  by  surface  pressure  measurements  increases 
gradually  along  the  entire  blade  surface  matching  very 
well  the  values  from  a  blade-to-blade  prediction 
program.  It  should  be  noted  that  the  local  flow  velocity 
becomea  supersonic  by  an  x/s  *  0.34.  The  predicted 
acceleration  parameter  drops  to  a  very  low  (favourable) 
value  very  early  on  the  blade  surface  and  becomes  only 
slightly  adverse  beyond  the  x/s  position  of  0,63.  This 
la  near  the  rearnost  point  where  wide  bandwidth 
transient  heat  transfer  data  are  collected  and  reported 
on  in  the  following  section.  It  is  also  near  the  point 
where  the  trailing  edge  shock  wave  from  the  adjoining 
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Figure  3  -  Baseline  Surface  Conditions. 

rotor  blade  is  expected  to  encounter  the  suctior 
surface.  Transitional  high  speed  data  is  not  reported 
past  this  point,  although  no  significant  change  in  the 
mean  heat  transfer  coefficient  was  observed.  The 
calculated  surface  curvature  is  also  shown  to  react  a 
low  value  a  short  distance  fron  the  leading  edge  and  tc 
remain  there  to  the  trailing  edge. 

natural  Transition  on  a  turbine  blade  surface  boundary 

LAYER 

Wide  bandwidth  heat  transfer  signals  from  surface 
thin  films  clearly  illustrate  the  important  difference 
between  the  iow  and  high  freestreac  turbulence  cases  as 
shown  in  the  transient  recorder  traces  of  Figure  4.  The 
heat  transfer  >ate  In  the  low  turbulence  case  (Run  572  3  ) 
remains  quiet  (laminar)  throughout  the  entire  measuring 
range  of  the  transient  data  (to  x/a  -  0.65).  The  high 
freestream  turbulence  case,  while  starting  eocewrat 
higher  than  the  laminar  case,  becomes  Increasingly 


Run  5713  Ti  <  0.8"  -  Flun  lu  s  IT 


Figure  4  -  Measured  Nusselt  Number.  The  Effect  of  Free- 
stream  Turbulence  on  Natural  Transition.  (The 
Vertical  Scale  for  Each  Trace  Range  from  0  to 
2000  in  Nusselt  Number.) 


dominated  by  sharp  transient  events  (consistent  with 
turbulent  spot  development,  amplification  and  gradual 
merger)  which  raise  the  heat  transfer  coefficient 
instantaneously  to  high  turbulent  levels  until  finally 
the  Nusselt  number  signals  become  core  characterised  by 
the  steady  turbulent  levels.  It  is  clearly  seen  at  the 
xfa  *  0.65  point  that  the  flow  is  at  the  turbulent  level 
more  than  half  of  the  time  dropping  precisely  to  the 
undisturbed  laminar  values  between  the  turbulent  events. 
This  process  starts  earlier  and  proceeds  faster  for  the 
higher  Reynolds  number  case  (not  shown  here). 


Figure  5  -  Natural  Transition  Heat  Transfer  Progression 
Along  the  Model  Surface. 


The  physical  process  of  turbulent  spot  development 
growth  and  merger  is  dramatically  illustrated  in 
Figure  5  in  which  a  time  synchronised  overlay  of  each 
instantaneous  heat  transfer  record  is  shown.  In  this 
figure  the  growth  and  rearward  convection  of  individual 
turbulent  spots  is  cleaily  seen  as  they  move  along  tbe 
surface,  again  with  the  levels  between  spots  dropping  to 
the  laminar  values.  This  growth  pattern  suggested  a 
procedure  for  the  quantitative  analysis  of  the  degree  of 
turbulent  activity  and  the  convection  rates  of 
individual  spots.  Figure  6  shows  expanded  time  records 
of  four  consecutive  positions.  By  selecting  a  turbulent 
threshold  level  slightly  above  the  laminar  value  the 
amount  of  time  that  the  digital  signal  is  above  this 
value  compared  to  the  total  time  of  the  trace  can  be 
used  to  define  an  interaittency  value  of  the  boundary 
layer  at  any  x/s  position  as 


n 


between  any  tim.  interval  t,  -  t,.  Also  clearly  seen  in 

Figure  6  is  the  regular  time  delay  for  the  turbulent 
spot  to  move  from  one  position  to  the  next.  If  one  were 
to  shift  the  signals  closer  together,  adjacent  heat 
transfer  records  would  be  seen  to  be  nearly  co-incident 
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Figure  6  -  Turbulent  Spot  Convection  Sate  and  Inter- 
oittency  Calculation. 

except  for  a  very  alight  but  regular  displacement  in  the 
direction  of  increasing  time.  Figure  7  presents  the 
value  of  the  intercii ttency  calculated  in  this  way  for 
the  two  conditions  of  design  and  1.5  x  design  Reynolds 
number  each  with  the  turbulence  grid  in  place.  The  rise 
of  the  interr.il ttency  above  zero  and  the  following  rapid 
rise  Is  seen  to  follow  for  some  time  the  predicted 
values,  reaching  a  value  near  unity  for  the  high 
Reynolds  number  case.  The  prediction  code  is  one 
developed  from  the  method  described  in  Patankar  and 
Spalding  (1970)  incorporating  a  model  for  mixing  in 
freestreao  turbulence  into  the  boundary  layer,  and 
allowances  for  curvature  and  changes  In  freestream 
conditions  during  transition,  outlined  by  Forest ( 1977 ) . 

The  convection  rate  for  turbulent  spots  was 
estimated  by  measuring  time  delays  for  individual  well- 
defined  spots  an'  from  cross-correlation  analysis  of 
adjoining  thin  film  signals.  Using  a  cross-correlation 
analysis 


Figure  7  -  Inters! ttency  Measurements  and  Predictions 
for  Natural  Transition. 


Figure  8  -  Example  of  Cross-Correlation  Analysis  o 
Adjacent  Heat  Transfer  Signals  in  a  Tran 
sitional  Boundary  Layer  (using  a  progr** 
developed  by  Doorly  (1983)). 


as  shown  in  Figure  8,  two  signals  (each  containing  102 
data  points)  are  compared  fcr  various  values  of  tit 
delay  and  the  peak  value  determined.  Correlatic 
coefficients  of  over  0.9  were  found  for  the  highl 
transitional  traces  late  on  the  blade  surface.  The  tim 
delay  wa«*  then  used  with  the  known  film  spacing  t 
estimate  a  mean  convection  velocity  of  the  discret 
events  associated  with  transition.  These  values  ar 
plotted  in  Figure  9  assuming  a  spot  origin  near  tt 
location  where  the  intermittency  rises  above  zero-  F> 
comparison,  the  spot  trajectories  are  compared  wi* 
curves  moving  at  various  fractions  of  the  bla, 
freestream  velocities.  The  mean  convection  rates  ar 
seen  to  be  between  0.7  and  0.8  of  the  local  freestrci 
velocity,  as  observed  in  low  speed  transitional  boundar 
layers . 
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Figure  9  -*  Measured  Convection  Rate  of  Turbulent  Spot 
by  Cross-Corre! at  1  on  Analysis  Cospared  wit 
Fractions  of  the  Local  Freestream  Velocity. 
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Ptgure  10  -  Comparison  of  the  Effect  of  Unsteady 
Wakes  and  Shock3  on  a  Freesfereaza  Hot-wire 
Output  with  Leading-Edge  Pressure  and  Heat 
Transfer  Rate  Measurements. 


NGV  WAKE  AND  SHOCK  WAVE  EFFECTS 

The  reaction  of  an  otherwise  laminar  boundary  layer 
to  a  higher  level  of  isotropic  freestream  turbulence  has 
been  detailed  above.  Clearly  transition  was  induced  by 
these  velocity  fluctuations  of  about  3%.  In  this 
section  results  are  presented  of  tests  at  the  design 
Reynolds  number  (0.92  x  10*  >  with  KGV  wake  and  shock 
structures  simulated  by  the  bar-passing  apparatus 
described  in  the  Introduction. 

The  intermittent  perturbations  to  the  inlet  flow 
caused  by  the  simulation  are  shown  in  the  high 
freestream  turbulence  case  in  Figure  10  in  terms  of  hot¬ 
wire  output  from  a  probe  mounted  in  the  freestream  inlet 
plane  at  mid-passage,  and  blade  surface  leading-edge 
measurements  of  pressure  and  heat  transfer  rate.  This 
data  is  normalised  with  respect  to  2  cycles  of  the  bar¬ 
passing  event,  so  that  it  is  possible  to  relate 
information  from  different  runs  in  terms  of  cycle 
fraction.  The  signals  show  a  background  turbulent  level 
extending  over  about  60%  of  the  cycle  characterised  by  a 
similar  type  of  signal  to  that  obtained  with  no  rotating 
bars*  All  three  signals  have  a  periodic  component  at 
bar-passing  frequency  with  two  characteristic  parts: 

<i)  Over  about  6%  of  the  cycle  rapid  changes  in 
level  of  the  order  of  5  to  10  ps  rise  and  fall  times  are 
observed.  The  pressure  signal  varies  by  25%  and  the 
heat  transfer  rate  by  approximately  ♦/-  50%.  This 
disturbance  is  attributed  to  shock  waves  generated  as 
the  tar  sweeps  pa3t  the  cascade.  By  examination  of 
Schlieren  photogr; .  hs  of  this  flow,  an  example  of  which 
is  given  as  Figure  11,  the  nature  of  the  NGV  simulated 
shock  structure  can  be  determined.  Clearly  two  shocks 
are  associated  with  one  bar-passing  event,  the  bow  and 
recotnpres3i on  shocks  that  wculd  be  expected  at  the  bar 
relative  Mach  number  in  steady  flow.  The  separation 
time  between  these  shocks  is  approximately  75ms,  marked 
as  in  Figure  10,  and  is  seen  to  ccrresponi  to  the 

time  interval  between  the  sharp  falls  In  2evel  to  the 
adjacent  sharp  peak. 

(il)  Following  these  rapidly  changing  events,  a 
second  less  marked  change  in  level  associated  with  the 
wake  is  also  visible  over  about  10%  of  the  cycle, 
denoted  33  the  "wake"  region  in  Figure  10, 


Figure  11  -  Schlieren  Photograph  Showing  the 
Intersection  of  the  Moving  Bar  Shock  Waves 
with  the  Turbine  Blade  Suction  Surfaces. 


It  is  possible  to  analyse  the  reaction  of  the  blad* 
boundary  layer  to  the  wake  and  shock  perturbations  by 
Investigation  of  the  sequence  of  tine-resolved  Musselt 
number  plots  feiven  in  Figure  12.  There  is  a  constant 
spacing  of  5  ma  between  the  thin  film  gauges,  starting 
with  gauge  3  at  an  x/s  value  of  0.12  through  to  gauge  15 
at  x/s=0.65  .  The  high  freestrean  turbulence  case  with 
wakes  and  shocks  present  (Figure  12(a))  is  markedly 
different  to  the  naturally  transitional  boundary  layer 
over  the  first  35%  of  the  surface,  with  similar  rapid 
rises  and  falls  in  Nusselt  number  to  the  perturbations 
evident  in  Figure  10.  This  shock  related  event  occurs 
on  the  early  suction  surface  due  to  a  shock /boundary 
layer  Interaction  starting  at  gauge  9,  close  to  the 
crown  of  the  auction  surface.  Examination  of  Schlieren 
photographs  (such  as  Figure  11)  indicates  that  the 

shocks  first  interact  with  the  boundary  layer  near  to 
gauge  9,  the  reflectl  n  point  coving  towards  the 
leading-edge  as  the  bar  moves  In  the  sare  direction. 
This  is  visible  on  the  early  gauges  on  Figure  12ia), 
occurring  first  on  gauge  9  then  moving  gradually  through 
gauges  7  and  5  and  finally  showing  on  gaut .  3.  The 
drop  in  surface  Nusselt  number  is  attributed  to  an 
unsteady  separation  and  the  rise  to  a  turbulent 
reattachment  both  caused  by  the  shock  boundary  layer 

Interaction.  The  effect  of  the  wake  is  not  clearly 

discernable  in  Figure  12(a)  and  to  aid  in  identification 
of  this  the  bars  were  rotated  at  a  lower  speed  such  that 
the  bar  relative  Mach  number  was  sutsrnie.  The  resul*? 


Pigure  12  -  Comparison  of  the  Effects  of  Shock  and  Wake- 
Passing  Eve  *  ‘  with  Undisturbed  Laminar 
Boundary  Layer  heat  Transfer. 


of  this  are  shown  as  Figure  12(b)  with  a  much  more 
clearly  identifiable  enhancement  in  heat  transfer  due  to 
this  wake.  This  extends  to  the  later  gauges  of  the 
surface  causing  the  boundary  layer  to  be  fully  turbulent 
over  the  extent  of  the  wake.  In  Figure  12(c)  the 
background  turbulence  was  reduced  to  less  than  0.8%  and 
the  periodic  disturbances  due  to  the  bar-passing  events 
are  more  clearly  evident.  The  early  suction  surface  has 
shock-related  phenomena  extending  well  into  the  cycle 
period  but  apparent  oscillations  in  Nusselt  number 
moving  with  the  shock.  Also  apparent  from  Figure  12(c) 
Is  the  intermittent  nature  of  the  turbulence  induced  in 
the  boundary  layer  by  the  wake  and  shock  interaction,  as 
the  boundary  layer  clearly  returns  to  its  undisturbed 
laminar  value  between  the  periodic  events.  The  heat 
transfer  enhancement  due  to  the  wake  :s  clearly  evident 
along  the  whole  surface. 

As  an  aid  to  quantifying  the  periodic  effects  due 
to  wake  and  shock  passing  phenomena  the  RMS  fluctuations 
caused  by  the  disturbances  in  the  coundary  layer  are 
plotted  for  each  gauge  in  Figure  13.  The  baseline  run 
(marked  as  x-x-x)  shows  a  constant  low  RMS  level  as 
would  be  expected,  and  the  higher  fre*»3tream  turbulent 
case  shows  a  steady  increase  following  the  boundary 
layer  intercl ttency  factor.  With  jus:  wakes  present  as 
was  shown  In  Figure  12(b)  a  general  enhancement  of  this 
RKS  value  is  evident  over  the  vhoie  surface,  whereas 
with  the  shock*  present  the  fluctuation  is  much  enhanced 
on  the  earlier  surface  but  returns  to  cr.dlsturbed  values 
towards  the  later  gauges.  The  values  come  close 
together  at  an  x/s  of  0.43  for  all  the  bar-passing 
cases.  Indicating  the  extent  of  the  shock  interaction 
regi on . 


it  *  -  V* 


Figure  13  -  Comparison  of  RMS  Nuaselt  Number  Flue 
tuations  over  the  Blade  Surface  for  Various 
Disturbances  at  the  Design  Reynolds  Number. 


CONCLUSIONS 

It  has  been  established  during  the  course  of  this 
study  that  the  Isentropic  light  Piston  Tunnel  facility 
combined  with  the  wide  bandwidth/high  sampling  rate  heat 
transfer  Instrumentation  has  proved  capable  of  tracking 
very  rapidly  progressing  unsteady  events  in  a  transonic 
boundary  layer.  Operating  under  a  simulated  unsteady 
gas  turbine  rotor  environment .  sensitive  detection  and 
precision  tracking  of  transient  shock,  wake  and  boundary 


layer  transitional  events  was  accomplished .  Detailed 
observation  of  the  heat  transfer  signals  confirmed  a 
classical  mod* 1  of  the  transition  process,  viz.  the 
sudden  appearance  of  sharp  spikes  of  turbulence  which 
multiply,  grow  and  track  downstream  at  a  rate  less  than 
the  freestream  velocity  until  gradually  merging  to  form 
a  completely  turbulent  boundary  layer.  The  tracking  of 
turbulent  spots  at  these  high  temperature  transonic  flow 
conditions  has  not  been  reported  before  to  the  best  of 
the  authors*  knowledge. 

A  second  major  finding  of  this  study  was  the 
observation  that  the  strong  unsteady  interaction  of  a 
double  shock  and  a  simulated  NGV  wake  with  the  rotor 
boundary  layer  did  not  have  any  measurable  "long  term" 
effects  apart  from  the  strong  excursion  in  heat  transfer 
associated  with  the  actual  passing  of  the  shocks  and 
wake.  The  heat  transfer  fluctuation  levels  were 
essentially  unchanged  far  removed  from  the  disturbance 
(in  time)  and  nearly  identical  at  the  rearmost  measuring 
point  except  for  a  turbulent  patch  associated  with  the 
shock,  wake  event  itaeir. 

The  interaction  of  the  shocks  and  wake  with  the 
rotor  establishes  in  more  detail  the  earlier  observation 
of  Ashworth  et  al.(1983)  and  Doorly  and  Oldfield  (1985a) 
of  strong  changes  in  local  heat  transfer  coefficient. 
The  tracking  of  the  interaction  over  the  surface  could 
be  followed  with  some  precision  with  the  time  resolution 
of  the  instrumentation  used.  Further  analysis  of  shock 
and  wake  interactions  Is  reported  by  the  authors  in 
Schultz  et  al.(I986)  along  with  some  prediction  models 
of  the  interaction  dynamics. 
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Tbs  strong  traiUng-edge  shook  waves  from  the  noxzle  guide  vena*  of  tranaonie  turbine  stags*  *an 
give  risa  to  interactions  with  tha  downstream  rotor  which  ara  significantly  aora  severe  than  la  tha  oasa 
with  lower  pressure  ratio  stages.  It  is  tharafora  important  to  study  suob  af facts  in  datall  both  fro* 
-b*  point  of  view  of  stags  power  output  and  mors  importantly  from  that  of  hast  transfer  ratas.  A  study 
h*s  baas  mads  of  a  transonlo  rotor  profile  in  a  static  cascade  in  which  tha  affect  of  shook  wave 
interaction  la  simulated  by  means  of  an  array  of  bars  rotating  at  tha  correct  spaed  and  spacing  upstream 
of  tha  stationary  rotor  blades.  Detailed  heat  transfer  rata  measurements  made  with  rapid  response  gauges 
enable  the  wake  and  shock  phenomena  to  be  separated. 


-  MCT  exit  velocity.  relative  bar  velocity 
*  Tangential  (or  true)  chord 
h  «  Thermal  conductivity  (V/mK) 

4  -  Heat  transfer  rate  Of/*1) 

H  -  isentropio  Mach  lumber  (based  on  local  statlo  pressure  and  inlat  total  praaaure) 

Ru  •  Rusaalt  Humber 

Re  *  Reynolds  Humber  (based  on  inlet  total  conditions,  isentropio  exit  Necb  lumber  end  tangential  chord) 
s  *  Blade  surface  perimeter 
t  *  Time 

T  -  Temperature  (C> 

Tu  -  Turbulence  level  -  u'/U 
u*  •  Fluctuating  velocity  (m/e) 
a  -  Telocity  (m/e) 

y  -  Rotor  relative  velocity,  casoede  inlet  velocity 

%  •  Surface  distance  from  the  leading  edge  stagnation  point 

jb  -  Distance  in  the  pitch-vise  direction  from  the  spanwlse  datum  position 

p  -  Cas  angle  (measured  from  the  axial  direction) 

Subscripts 
«  -  Freestresm 

o  *  Total 

1.2  »  Inlet.  Outlet 

mess  *  Measuring  point 

h.m.t  ■  Hub.  mean,  tip  (or  tangential  in  chord  definition) 
rel  *  Relative  bar  condition 
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edge.  The  ouotdi  baa  a  spaa  or  50  mi  at  blade  lolat  and  56.095  aa  at  exit  with  the  expansion  on  one 
side  only  aa  illustrated  In  Pig.  t(a).  A  turbulence  grid  20S  am  upstrean  of  tbe  cascade  provides  a  level 
of  u’/O  of  approximately  3%.  Ope  tree*  and  downs  trees  static  and  total  pressures  were  seaaured  routinely 
to  establish  the  correct  cascade  operating  conditions  and  are  reported  in  sore  detail  by  Ashworth  et 
el*1*.  Tbe  operating  conditions  st  the  nominal  engine  design  point  are  given  below. 

The  heat  tranafer  gauges  used  In  this  study  are  oonventlonal  thin  filn  surfaoe  resistanoe 
thermometers  widely  used  for  tbe  deterelnatloo  of  beet  transfer  rate  in  short  duration  facilities1*. 
Data  fros  22  auoh  gauges  were  stored  In  a  digital  transient  recorder  seapling  up  to  16  channels  st 
SOD  kHx  for  esob  ohannel  or  were  input  dlreotly  to  tbe  A/D  oonverter  at  a  slower  rate  of  400  Hz  for  sose 
of  the  64  available  A/D  channels  when  tlse  average  date  only  were  required,  also  used  for  sessuresenta 
such  ss  Inlet  end  exit  statlo  pressures.  The  looetions  of  these  best  transfer  gauges  on  tbe  blade  are 
shown  in  Pig.  2(a)  and  given  in  terns  of  the  surface  length  ’x*  to  periaeter  ♦s’  from  the  stagnation 
point.  A  nore  detailed  study  of  the  reaction  of  the  suction  surface  boundary  layer  to  both  freest ream 
turbulence  and  tbe  wake-passing  phenomena  is  also  reported  below.  For  this  study  another  blade  was 
instrumented  with  thin  film  gauges  which  were  only  4  na  long  aa  compared  to  10  as  for  tbe  previous  tests, 
and  were  nore  olosely  epeced  around  tbe  profile,  as  shown  in  Pig.  2(b).  The  nodel  points  in  Pig.  2(b)  are 
those  referred  to  in  Pig.  6.  Surface  pressures  were  aeesured  in  e  tine-averaged  manner  using  Sensya 
semiconductor  transducers  type  LX-1620D  operating  In  e  differential  node  as  reported  in  nore  detail  by 
Ashworth  et  el.1*.  Baseline  experimental  results  have  been  reported  by  Ashworth  et  al.  but  for  clarity 
acne  of  this  data  is  referred  to  in  the  present  paper.  The  velocity  triangle  for  tbe  steady  nominal 
design  case  is  shown  in  Fig.  Kb)  which  includes  the  effect  of  the  reduced  NGV  exit  wake  velocity  C£  on 
the  rotor  inlet  angle 


HMtt  Rmm’  Imiflt  Without  Ml  Infraction 


A  comparison  of  baseline  data  with  no  rotor/wake  interaction  is  given  in  Pig.  3(a)  for  tbe  two  cases 
of  low  (<  0.1%)  freest  ream  turbulence  and  with  a  turbulence  level  of  spproxinately  3%.  It  will  be  seen 
that  tbe  turbulence  generated  by  tbe  bar  grid  le  sufficient  to  bring  the  region  of  boundary  layer 
transition  forward  fros  about  60%  x/s  on  the  pressure  surface  and  50%  x/s  on  the  suction  surface  to  10% 
and  20%  respeotively .  All  of  this  data  was  taken  at  the  nominal  design  cascade  operating  conditions; 


432  K 

H  ■  i.l* 
exit 


total 


(based  on  Isentroplc  Inlet  total  and  pltcbwise  averaged  exit  static  pressures) 

He  -  0.919  x  10* 

-  51.06*. 

Tbe  circled  numbers  refer  to  heat  tranafer  gauges  identified  in  Pig.  21a).  The  heat  transfer  rate 
is  presented  In  terms  of  a  non-dimensional  Nusselt  Humber.  defined  as: 
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Instantaneous  and  Mean  Heat  Transfer  with  Wake  Interaction 

The  heat  transfer  rate  to  the  blade  with  the  additional  effect  of  wake  interaction  is  illustrated  in 
Figs.  3(b)  and  3(c)  for  both  cases  of  effectively  zero  freeatreaa  urbulence  and  -  3%  u’/U.  Prom 
Pig.  3(c)  it  will  be  seen  that  there  le  an  overall  iocraaaa  in  heat  transfer  rate  over  both  the  pressure 
and  auction  surfaces.  The  pressure  surface  heat  transfer  is  enhanced  over  practically  the  whole  length 
although  the  dominant  effect  Is  observed  for  values  of  x/s  <  70%.  On  the  pressure  surface  the  effect  of 
wake  interaction  persists  to  about  x/s  *  30%.  Examples  of  instantaneous  heat  transfer  rates  are  inset  In 
the  Figure  end  e  more  extemslve  'atlas'  of  results  is  >tven  lb  Ashvortb  et  el.1*.  The  heat  tranafer 
rates  with  waxe  interaction  and  with  effectively  zero  freeetreas  turbulence  are  shown  in  Pig.  3(b).  As 
expected  there  is  e  marked  increase  in  the  level  of  heat  transfer  rate  over  almost  tbe  entire  pressure 
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(a)  natural  Transition  of  the  Suction  Surface  Boundary  Layer 

Vide  bandwidth  heat  transfer  signals  obtained  from  the  surface  thin  file  gauges  clearly  illustrate 
the  important  differences  between  the  low  and  high  freestreas  turbulence  oases  ss  above  In  the  high 
frequency  traces  of  Fig.  6(a).  There  is  a  constant  spacing  of  5  an  between  the  thin  fiia  gauge  results 
shown,  starting  with  gauge  3  at  an  x/s  value  of  0.12  through  to  gauge  15  et  x/a  »  0.65.  The  low 
turbulence  oase  (Run  3723)  remains  quiet  (laminar)  throughout  the  entire  measuring  range  of  the  transient 
data  <to  x/s  -  0.65).  The  surface  best  transfar  for  tbs  high  freestrsaa  turbulence  case,  while  starting 
somewhat  higher  than  the  laminar  case,  becomes  increasingly  dominated  by  sharp  transient  events 
(consistent  with  the  theory  of  turbulent  spot  development,  growth,  end  gredual  aerger  as  proposed  by 
Eatons**  and  verified  by  aany  others  (e.g.  Schubauer  and  Klabanoff**) .  These  spots,  which  raise  the  heat 
transfer  coefficient  instantaneously  to  high  turbulent  levels,  continue  to  grow  and  merge  until  finally 
the  Nusselt  Nuaber  signals  becoae  increasingly  characterised  by  the  "steady”  turbulent  levels.  It  la 
clearly  seen  that  by  the  x/s  -  0.65  station  the  flow  is  at  the  turbulent  level  more  than  half  of  the  time 
but  drops  precisely  to  the  undisturbed  laainsr  values  between  the  turbulent  events.  In  other  tests 
conducted  st  1.5  x  Be  design,  the  boundary  layer  was  fully  turbulent  at  this  looatlon. 

The  physical  process  of  turbulent  spot  breakdown  to  turbulance  can  be  seen  In  Fig.  6(a).  The  growth 
and  rearward  convection  of  Individual  turbulent  spots  is  clearly  seen  as  they  nova  along  tba  blade 
surface.  The  spot  signals  grow  in  height  and  width  as  the  spots  cover  aore  of  each  succeeding  thin  file 
gauge  and  shift  in  time  downstream.  This  breakdown  process  Is  quantified  In  aore  detail  in  Ashworth** 
where  interalttenoy  levels  are  estimated  froa  the  digital  time  records  and  spot  oonvectlon  rates  are 
estimated  from  cross-correlation  analysis  of  adjoining  thin  film  signals. 

(b)  Detailed  Nake  and  Shook  Interaction  Effects 

It  is  possible  to  analyse  the  reaction  of  the  blade  boundary  layer  to  the  wake  and  ahook 
perturbations  with  2  bars  rotating  by  investigation  of  the  sequence  of  time-resolved  Nusselt  Nuaber  plots 
given  In  Fig.  6(b).  (o)  and  (d) .  The  high  freeatreaa  turbulence  case  with  wakes  and  shooks  present 
(Fig.  6(b))  is  markedly  different  to  the  naturally  transitional  boundary  layer  (Fig.  6(a))  over  the  first 
35%  of  the  surface,  with  similar  rapid  rises  and  falls  in  Nusselt  Nuaber  to  the  perturbations  evident  in 
Fig.  3.  This  shook  related  event  occurs  on  the  early  suction  surface  due  to  a  shock/ boundary  layer 
interaction  starting  at  gauge  9,  olose  to  the  crown  of  the  suction  surface.  Examination  of  the  Sohlleren 
photographs  (Fig.  7)  indicates  that  the  shooks  first  Interact  with  the  boundary  layer  near  to  gauge  9. 
the  reflection  point  moving  towards  the  leading-edge  as  the  bar  aoves  in  tbe  saae  direction.  This  is 
visible  on  the  early  gauges  on  Fig.  6(b),  occurring  first  on  gauge  9  then  moving  gradually  through  gauges 
7  and  5  and  finally  showing  on  gauge  3.  The  rapid  drop  In  surface  Nusselt  Nuaber  is  attributed  to  an 
unsteady  separation  and  the  rise  to  a  turbulent  re-attachment  both  caused  by  the  shock  boundary  layer 
interaction.  Tbe  effect  of  tbe  wake  is  not  clearly  dlsQernable  in  Fig.  6(b)  and  to  aid  in  identification 
of  this  the  bars  were  rotated  at  a  lower  speed  such  that  the  bar  relative  Mach  Nuaber  wee  subsonic.  The 
results  of  this  ere  shown  as  Fig.  6(c)  with  a  ouch  aore  clearly  identifiable  enhancement  in  heat  transfer 
due  to  this  wake.  This  extends  to  the  later  gauges  of  the  surface  causing  the  boundary  layer  to  be  fully 
turbulent  over  the  extent  of  tbe  wake.  In  Fig.  6(d)  the  background  turbulence  was  reduced  to  less  tban 
0.6%  and  tbe  periodic  disturbances  due  to  the  bar-passing  events  are  aore  clearly  evident.  The  early 
auction  surface  has  shock-related  phenomena  extending  well  into  tbe  cycle  period  with  apparent 
oscillations  In  Nusselt  Number  moving  with  the  shock.  Also  apparent  from  Fig.  6(d)  Is  the  intermittent 
nature  of  tbe  turbulence  induced  in  the  boundary  layer  by  the  wake  and  shock  interaction,  as  the  boundary 
layer  olearly  returns  to  its  undisturbed  laminar  value  between  the  periodic  events.  The  heat  transfer 
enhanoeaect  due  to  the  wake  is  clearly  evident  along  the  whole  surface. 

In  summary.  It  appears  that  the  state  of  the  turbine  boundary  layer  seems  to  bs  controlled  by  the 
level  of  freeatreaa  turbulence  except  in  the  region  where  tbe  actual  shock  and  wake  peases  through  tbe 
cascade  paasage. 


(i)  A  prediction  of  the  flowfield  velocity  is  Bade.  la  this  esse  using  the  Denton  scheme*4  end  stored 
by  the  program. 

(11)  Op^ia  calculated  (assumed  constant  across  the  span),  the  center-line  of  the  uadi* tor ted  wake  is 
calculated  froa  the  specified  bar  position  and  the  width  added. 

<ili)  The  wake  is  shifted  back  In  time  so  that  tha  bar  will  return  to  its  correct  position  following  the 
■arobing  process  of  the  prediction  routine. 

(lv)  Prom  this  initial  position  elements  of  the  wake  are  contacted  by  small  tlae  steps  using  the  local 
velocity  interpolated  fro*  the  prediction  until  the  bar  reaches  tha  specif  lad  location.  The 
differential  velocities  in  the  flowfleld  cause  distortion  of  the  wake  along  its  length  and  across 
its  width  as  it  is  accelerated  through  the  passage. 

This  staple  scheme,  which  could  be  incorporated  in  the  blade  design  process,  agrees  well  with  the 
positions  of  the  wake  on  the  Scblieren  photographs,  so  that  the  fraction  of  tiae  that  tha  heat  transfer 
rate  to  the  surface  is  affected  by  the  wake  could  be  calculated  and  included  in  the  interaittency  tern  of 
a  prediction.  It  also  demonstrates  that  to  the  level  of  our  aeasuring  ability  sacond-ordar  effects  such 
as  the  "negative  Jet"  effect  do  not  significantly  affect  the  wake  position,  aa  referrad  to  in  Doorly*1, 
and  that  the  assumption  that  this  flow  unsttadinesa  may  be  superimposed  on  the  ateady  flowfleld  is  m 
valid  approximation. 

Stock  IrtfllsUaa 

Although  it  was  possible  to  aodel  the  wake  in  2-dimensional  terms,  the  shock  structure  associated 
with  the  transonic  nature  of  the  bar  (as  Mp#i  varies  from  1.06  et  the  hub  to  1.25  at  the  tip)  Is  not 
strictly  2-dimensional,  as  has  been  seen  in  the  Schlieren  photographs  (Fig.  7).  An  attempt  was  aade  to 
predict  the  positions  of  the  bow  and  recompresaion  shocks  due  to  the  bar  in  order  to  allow  trajectory 
rate  calculations  to  be  aade  relating  to  the  heat  transfer  measurements,  and  to  aid  in  understanding 
phenomena  observed  in  the  Schlieren  photographs.  An  example  of  such  a  prediction  is  presented  as  fig*  9. 
The  prediction  is  based  on  the  assumption  that  for  a  small  element  of  the  ber  the  flow  is  2-dlmensional 
with  respect  to  the  plane  containing  0  ^  and  normal  to  tha  bar  axis,  a  quasi -3 -dimension*!  *ppro*cb 

described  In  detail  in  Ashworth*1.  As  U^^varles  in  both  magnitude  and  direction,  so  doea  the  associated 
shock  structure,  and  unlike  the  wake,  this  varietlon  should  he  accounted  for.  The  method  of  prediction 
was  computerised  as  follows: 

(i)  For  assumed  constant  inlet  conditions  and  a  specified  position  of  the  bar  Mp<1  is  oalculated  and 
its  direction  determined. 

(11  >  The  equation  of  the  bow  shock  in  the  2-dimenslonal  plane  described  above  is  derived  from  M  ^  using 
the  method  described  in  Shapiro**,  shook  properties  and  turning  angles  obtained  by  curve-fit 
equations  where  necessary  over  a  range  of  Mach  Number  from  1.0  to  1.5. 

(iii)  The  recompression  shock  is  assumed  to  be  straight  and  inclined  et  the  Mach  angle  (sin"  1/Mre,>  to 
the  direction  of  M  ^ ,  with  a  virtual  origin  two  diameters  downstream  of  the  bar. 

(iv)  The  shocks  are  chopped  and  refracted  if  they  are  downstream  of  the  blade  axial  leading-edge  point, 
with  allowance  made  for  regeneration  of  the  shook  as  it  moves  away  from  this  point. 

<v)  Finally  the  points  of  intersection  with  the  blade  suction  surface  are  calculated,  and  tha  shocks 
are  simply  reflected  from  an  origin  mid-way  between  the  two  intersection  points. 

Comparison  with  the  Schlieren  photographs  is  encouraging,  the  position  shown  in  Fig.  9  corresponding 
to  time  2  In  Fig.  7.,  and  it  is  hoped  that  information  from  this  simplified  model  will  prove  useful  as  an 
aid  to  understanding  the  complex  shock  movements  in  turbine  passages,  in  spite  of  the  many  simplifying 
assumptions  aade  In  this  prediction  (such  as  allowing  for  no  variation  of  the  freeatream  velocity). 

CcMltnUtti 

It  baa  been  es  abliahed  during  the  course  of  this  study  that  the  leentropio  tight  Piston  Tunnel 
facility  combined  with  the  wide  tandwidth/high  sampling  rate  heat  transfer  instrumentation  has  proved 
capable  of  tracking  very  rapidly  progressing  unsteady  events  in  a  transonic  boundary  layer.  Operating 
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Fig.  1(a)  Arrmgcwnt  of  the  rotating  bar  vakc 
generator  and  caecede. 


Fig,  l(b)  Velocity  trlanglea  at  inlet  for  the 
steady  nominal  daalgn  case  at  mid* 
span,  shoving  the  reduced  incidence 
caused  by  the  deficit  In  the  vake. 
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Fig*  2(a)  Co-ordinates  of  the  original  heat  traoafar  gauges  on  the  blads  profile. 

Fig*  2(b)  Co~ordinates  of  the  heat  transfer  gauges  for  the  detailed  suction  surface  study, 
referred  to  la  Fig.  6. 
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Fl| .  3(a)  Effect  of  fraaatraaa  turbulence  level  on  Been  beet  tranefer  race  without  vakee  for 
the  noainal  daalgn  caae.  “  1.18,  Re  ■  0.919E6. 
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3(b)  Effect  of  wake  and  shock  lntsractton  on  assn  heat  transfer  rate  at  the  high  free- 
atreaa  turbulence  case  ( “  IX)  and  design  operating  conditions.  Inset  are  typical 
transient  recorder  algnala  for  the  26  bar  experiment  at  two  x/e  locations  shoving 
e  1  as  Interval  of  the  record  together  vlth  the  tlae-averaged  values  without  the 
wake  and  shock  Interactions. 
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Pt«.  3(c)  Effect  of  »e he  and  ehock  Interaction  on  naan  haat  tranafar  rata  for  tho  low  frao- 
atraaa  turbulaaca  caaa  (<  .81).  othenrlaa  aa  U  ft».  3(b), 
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§1  RUN  CONO  SYMBOL  I  GRID  WAKES 
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suction  surface  fractional  perimeter  pressure  surface 
HE  L/E  T/e 

Fig.  Me)  Effect  of  Incidence  variation  on  mean  beat  transfer  rate  et  the  noBlnel  design  case 
(»2  *  i . 18,  Re  •  0.919E6)  with  no  wake  and  shock  Interaction. 


Fig.  Mb)  Effect  of  Incidence  variation  on  mean  surface  isentropic  Mach  Number  at  the  nominal 
design  caee  (M2  ■  i.lfc.  Re  *  0.9L9E6). 


CASCADE  SIMULATION 


ROTOR  FI1ED 

mi  mqv  1  mg 


rig.  5(a)  Illustration  of  the  cooperative 

velocity  triangles  for  the  engine 
and  cascade  simulation  of  the  MCV/ 
rotor  interaction. 


Fig.  5(b)  He stu resents  of  tha  unataady  distur¬ 
bance  it  inlet  by  a  hot-wire  in  the 
freeatreaa  and  hast  transfer  and 
pressure  measurements  at  the  stag¬ 
nation  point. 


Fig.  6  Detailed  time-resolved  heat  transfer  rate  measurements  to  the  blade  suction  surface.  The  dotted  lines  are  the  baseline 
undisturbed  results  for  the  low  freestream  turbulence  case  (<  IX,  Run  5723). 
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Fig,  9  Predicted  ehock  position  corresponding  to  eiaa  2  of  Fig.  7 .  with  th*  bar  Inclined  et 
*  2*  to  Ite  datua  position.  Only  the  helf  of  the  ehock*  dovnetreea  of  the  bar  are 
shown  for  clarity.  The  rtfrsctlon  and  reflection  of  both  ehocke  can  be  seen,  with  the 
re-coapreaslo«  ehock  aore  2-dlaenslonaI  in  fora  then  the  detached  bov  shock*  due  to 
•t end-off  distance  variation  with  bar  relative  Mach  Nuaber.  Both  shock*  era  weaker 
tower da  the  hub  due  to  the  lover  Mach  Nuabtr  there,  with  the  bov  ehock  strength  leas 
than  the  re-covpreeslon  ehock  strength  at  each  radius. 
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